Nanowire metamaterials with extreme optical anisotropy 
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We study perspectives of nanowire metamaterials for negative-refraction waveguides, high- 
performance polarizers, and polarization-sensitive biosensors. We demonstrate that the behavior 
of these composites is strongly influenced by the concentration, distribution, and geometry of the 
nanowires, derive an analytical description of electromagnetism in anisotropic nanowire-based meta- 
materials, and explore the limitations of our approach via three-dimensional numerical simulations. 
Finally, we illustrate the developed approach on the examples of nanowire-based high energy-density 
waveguides and non-magnetic negative index imaging systems with far-field resolution of one-sixth 
of vacuum wavelength. 



The anisotropy of effective dielectric permittivity is 
widely used in optical, infrared (IR), THz and GHz 
sensing, spectroscopy, and microscopjsii^i^ii^. Strongly 
anisotropic optical materials can be utilized in non- 
magnetic, non-resonant optical media with negative in- 
dex of refraction, and have the potential to perform 
subdiffraction imaging and to compress the radiation 
to subwavelength areas^i^i^ii. The performance of these 
polarization-sensitive applications can be related to the 
relative difference of the dielectric constant along the dif- 
ferent directions. In the majority of natural anisotropic 
crystals this parameter is below 30%^. While it may 
be sufficient for some applications, a number of excit- 
ing phenomena ranging from high-performance polariza- 
tion control'' to subwavelength light guidingSii^i^ to planar 
imaging^ require different components of a permittivity 
tensor to be of different signs. 

In this Letter we study the perspectives of using 
nanowire composites as meta-materials with extreme op- 
tical anisotropy. We demonstrate that even 10% stretch- 
ing/compression of the nanowire structures may dramat- 
ically affect the electromagnetic properties of these sys- 
tems and change the sign of components of the permittiv- 
ity tensor. We present an analytical description of wave 
propagation in anisotropic nanowire composites - Gener- 
alized Maxwell-Garnett approach (GMG), and verify our 
technique via three-dinrensional (3D) numerical simula- 
tions. Finally, we illustrate our approach on the exam- 
ples of several nanowire-based systems for light compres- 
sion below the diffraction limit and negative refraction- 
imaging with far-field resolution of Ao/6 (with Aq being 
free-space wavelength). 

The use of metallic wire mesh as anisotropic low- 
frequency plasma has been proposed in^ and experimen- 
tally realized for normal light incidence in-'--. However, 
the applicability of these nanowire-based materials for 
any non-trivial geometry involving oblique light incidence 
or wave-guiding is still considered to be questionable due 
to strong nonlocal intcractions'-i, that may potentially 
result in positive components of the permittivity tensor. 
Furthermore, the majority of existing effective-medium 



theories (EMTs) '^'^^i-'^'-''^ are limited to the optical re- 
sponse of nanowires that are isotropically distributed in 
the host material. The predicted response of these sys- 
tems is almost independent of nanowire distribution and 
is described by a single parameter - nanowire concen- 
tration. These existing techniques are therefore not ap- 
plicable for practical composites where the geometry is 
anisotropic due to fabrication process or as a result of 
a controlled mechanical deformation^^. Understanding 
the optical behavior of nanowire structures beyond one- 
parameter EMT is the main purpose of this Letter. 

The geometry of the nanowire composites considered 
in this work is shown in Fig[TJ The nanowires with per- 
mittivity e'" are embedded into a host material with 
permittivity e°"*. The nanowires are aligned along the 
z direction of Cartesian coordinate system. We as- 
sume that the nanowires have elliptic crossections with 
the semi-axes Tx and ry directed along x and y coor- 
dinate axes respectively. We further assume that the 
homogeneous nanowire conrposite may be compressed or 
stretched, leading to the anisotropic distribution of indi- 
vidual nanowires. The typical separations between the 
nanowires x and y directions are denoted by Ix and Zjp^. 
In this work we focus on the case of "homogeneous meta- 
material" , when the inhomogeneity scale is smaller than 
the wavelength {Tenia ^ Aq) and nanowire concentra- 
tion N is small so that the lattice feedback effects can be 
treated using perturbative techniques as described below 
[for the case of e'" < additional requirenrent ra < cr, 
with a being skin depth in wires, must be fulfilled]. Un- 
der these conditions. Maxwell equations have free-space- 
like solutions that can be represented as a series of plane 
electromagnetic waves propagating inside a material with 
some effective dielectric permittivitj*ii e : 
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The angular brackets in Eq.(IT]) denote the average over 
microscopically large (multi-wire), macroscopically small 
(subwavelength) region of the space, with Greek indices 
corresponding to Cartesian components, and assumed 
summation over the repeated indices. If both e'" and e™* 




the effective permittivity e°*^, e°^ and e°f . Using the con- 
tinuity of the Ez component, Eq.(IT]) yields: 
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Note that similar to what was found in Refsiii'i^iii^'i^, the 
single parameter that determines the z component of the 
permittivity in the effective medium regime is nanowire 
concentration N. 

To find {x, y) components of the e*''^ we use the 
Maxwell-Garnett (MG) technique^ii^iii. This approach 
assumes A^ ^ 1 so that the local field in the composite is 
homogeneous across a nanowire. The fields D and E are 
then averaged over a typical nanowire cell, and Eq.(IT]) is 
used to extract the effective permittivity of a material. 
Naturally, average fields will have two contributions: one 
coming from the fields inside nanowires E™, and the sec- 
ond one coming from the fields between nanowires i5°"*. 
The derivation of an EMT is therefore equivalent to un- 
derstanding the relationship between E™,E°^^ and the 
external field acting on the system E^ . Conventional MG 
approach assumes that E""^^ — iJ°— iH, which is true only 
for the case of isotropic nanowire distributions. The cru- 
FIG. 1: (color online) (a) Schematic geometry of a nanowire ^ial point of this work is that both E'"" and £;°"* are 




composite, (b) e^" for e'" = -2.5, e°" = 2 (Ag nanowires 
in a polymer for Ao — SBOnin). (c-e) eJrj, (red triangles, solid 
lines) and ej^* (orange rectangles, dashed lines) for the com- 
posite in (a), with K^ ^ Q,o, = Q (c), Kx = 0.2; 0.^ = -0.2 (d). 
Ax = 0.2; ilx ~ (e), A^ = 0;Qx ~ —0.2 (f); quasistatic nu- 
merical calculations (symbols); and GMG (lines); dotted lines 
in (d)- ■ ■ (f) are identical to lines in (c). Bottom insets show 
crossections of composites for A'^ = 0.35; Black dashed lines 
in top insets illustrate the effect of small losses in nanowires 
(e'" = —2.5 -|- 0.04J); note that this effect is negligible for 
l^effj < 20. The breakdown of GMG occurs when the local 
field becomes inhomogeneous on the scale of Va 



are isotropic, the normal axes of the tensor of effective 
dielectric permittivity will coincide with the coordinate 
axes. Thus, in the selected geometry the permittivity 

Sapefp, with Sap being 



eff 

a(3 



tensor becomes diagonal: e 
the Kronecker delta function. 

We now derive the expressions for the components of 



strongly influenced by the nanowire distribution (given 
by the parameters lx,ly), and nanowire shape (described 
by rx,ry) along with nanowire concentration TV. 

We now derive the analytical expressions for E™ 
and E°^'^ . The typical excitation field acting on a 
nanowire in the sample will contain the major contri- 
bution from external field E'^ and the feedback field 
scattered by all other nanowires in the structure xE'^ , 
resulting in the effective excitation E^ + ^ x' E^ = 
[5ap — Xap\''^E^- For rectangular, triangular, and other 
highly-symmetrical lattices, as well as for a wide-class of 
random nanowire micro-arrangements, the feedback ten- 
sor becomes diagonali^, with the effective field acting on 
a nanowire being [1 — Xaa\~^E'^^. 

Using the dimensionless function S{^) — ^^ ■ ^ij^^fi -m 
with summation going over all pairs of i, j except coor- 
dinate origin, the summation of 2D dipole fields over 
rectangular lattice shown in Fig[TI^a) yields^ 
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where we introduced the lattice distortion vector 
{Ax, Ay} — {Ix/ly — ^,lyl^x ~ 1}, and polarization term 



na(e" 



'^)], with {ux, Uy} = {ryl(rx + 
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Note that the feedback parameter vanishes only for 



,)\ being the depolarization factor a^'*'^^ . 



isotropic nanowire distribution Ix — ly, corresponding 
to the well-known MG result^'ii'i^ii^'i^. 

This inter-wire interaction changes the "microscopic" 
field acting on the individual nanowires, and thus it 
directly affects both (homogeneous) field inside the 



nanowire E™—, 
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and the average field in-between the nanowires i5°"* . Di- 
rect calculation of the average dipole^^ E°^^ of a given 
inclusion over the typical meta-material cell yields: 
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with Q{N) — TT ~ 1 — N{tt — 2) and shape vector 

Combining Eqs. (|ll4l5p we arrive to the following ex- 
pression for the in-plane components of permittivity in 
GMG approach: 
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FIG. 2: (color online) (a) dispersion of the fundamental TM 
(red triangles, solid lines) and TE (orange squares, dashed 
lines) modes in a waveguide with d = 400nm with nanowire 
composite core; e'" = 13; e°"* = -120; A^ = 0.2; 0^ = 
—0.2; Ix = 40nm; r^ = lOnm; Note the negative refraction 
mode, predicted ir^; (b) modal propagation constant for 
Ao = 1.5/^771 as a function of d; numerical solutions of 3D 
Maxwell equations (symbols); Eas. (|2l6l7p (lines). The break- 
down of GMG correspond to I/cq/q,! ~ 1 

To verify the accuracy of the developed GMG tech- 
nique, we generate a set of nanowire composites with 
given values of N, e'", e°"', A, and fl, excite each com- 
posite with a homogeneous field, and use the commercial 
finite-element partial differential equations solver, COM- 
SOL Multiphysics 3.2^^ to solve Maxwell equations, find 
the microscopic filed distribution, and calculate the aver- 
age values of E, and D over the volume of a composite, 
yielding e . In the simulations we used both random and 
periodic nanowire composites; the number of nanowires 
was sufficient (> 10^) to eliminate the dependence of 
e"*^ on the sample size [finite-sample-size artifacts] . FigU] 
shows the excellent agreement between GMG approach 
presented in this work and numerical solution of Maxwell 
equations in quasistatic limit for concentrations N < 0.3 
and deformations 117^1, |Aq| < 0.3. Our simulations in- 
dicate that the quasi-static material properties are fully 
described by average parameters {N,la,ra). This par- 



ticular property of the effective-medium composites indi- 
cates high tolerance of anisotropic mctamaterials to pos- 
sible fabrication defects. 

As expected, the field distribution across the nanowire 
structure and e^^ are strongly affected by TV, as well as 
A and fl. Thus, even 10% anisotropy in inclusion shape 
or distribution may result in change of sign of dielectric 
permittivity. Such an effect opens the possibility to create 
optical materials with widely controlled opto-mechanical 
properties, potentially leading to new classes of polar- 
izers, beam shapers, polarization-sensitive sensing and 
fluorescence studies, as well as for a wide class of ultra- 
compact waveguides^"*' since the material properties may 
be tuned between e « and |e| 3> 1. Some of these 
applications are described below. 

As it has been noted for GHz systems in^ -, the com- 
ponents of e°^ may be strongly affected by the spatial 
dispersion. To clarify these effects we used COMSOL 
package to identify the eigen waves propagating in x- 
direction through a planar waveguide with a compos- 
ite core consisting of a rectangular array of 10% Ag 
nanowires in Si host, extending from z — to z — d (see 
Fig. la), bounded by perfectly conducting walls (see-^ for 
the detailed explanation of the effects of waveguide walls 
and material absorption on the mode propagation). In 
FigEfa) we show the agreement of the results of numeri- 
cal solutions of 3D wave equations with the EMT dynam- 
ics of TE and TM modes propagating in a waveguide with 
homogeneous anisotropic core, described by 
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with oj — 2tt/Xo and kx,c,d being the modal wavevector, 
speed of light in the vacuum, and waveguide thickness, re- 
spectively. Note that this system does not support TEM 
modesi. 

FigHUb) illustrates one of the applications of nanowire- 
based optical composites, high-energy-density waveguide 
- a subwavelength structure supporting propagating vol- 
ume modes. It is important to point out that in contrast 
to uniaxial media, anisotropic (e®*^ ^ 6°^) nanowire com- 
posites can simultaneously support both n > 0-TE and 
n < 0-TM waves {n = kxc/oj). Moreover, the in-plane 
anisotropy (induced, for example, by deformation) can be 
used as a controlling mechanism in nanoscale nanowire- 
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FIG. 3: (color online) Planar waveguide imaging via nanowire 
materials Aq — 1.5/im; d — 360nm; n > material: e = 13; 
nanowire composite (4 < 2 < 8): r^ ~ Vy = lOnm; Ix = ly ~ 
50nm; e'" = —120; e°"* — 2. (a) intensity across the system; 
(b) intensity distribution across the focal plane; insets show 
planar lens geometry and dispersion of a negative-index mode; 
dash-dotted line corresponds to d = 360nm (a). 



based pulse-management devices. 

It is clearly seen that the propagation of these modes is 
adequately described by GMG technique when |A:q^q| ^ 
1. As expected, the material properties in EMT regime 
are independent of nanowire micro-arrangements (type of 
crystalline lattice), while the exact point of EMT breakup 
(|fca^a|max) depends on local geometry and is maximized 
for almost-rectangular lattices [assumed in derivation of 
Eq.([3])]. Note that the real requirement for EMT appli- 
cability, IfcctZal <^ 1, is different from the commonly used 
criterion la <C Aq. Indeed, our simulations show that 
the spatial dispersion leads to cut-off of the modes with 
l^a^al ^ Ij similar to what has been predicted for GHz 



wire systems -''^ and nanolayer-based photonic funnels^. 

Another application of nanowire structures, non- 
magnetic negative-index materials^ is illustrated in FigO 
It is seen that the nanowire materials may be used to 
achieve sub-diffraction (Ao/6) far-field resolution in the 
planar-lens geometry. 

In conclusion, we have developed the effective-medium 
theory (GMG) that adequately describes the opti- 
cal properties of nanowire composites with anisotropic 
crossections and arrangements. Limitations of the pro- 
posed approach have been studied via numerical mod- 
eling. We demonstrated that the nanowire composites 
can be used to achieve extreme anisotropy at optical 
and IR frequencies, with controlled effective permittivity 
ranging from e^— ltoe«Otoe^l- thus lead- 
ing to practical implementations of high-energy-density 
waveguides^, novel polarization-sensitive detectors, and 
recently proposed non- magnetic negative index systems^. 
Finally, we note that the technique presented here can 
be readily applied to dielectric, plasmonic, and polar- 
wire composites at optical, IR, and THz frequencies, and 
can be further extended to the cases of non-aligned inclu- 
sions, anisotropic e'" and e°"', and 3D composites similar 
to what have been done for isotropic-arrangement cases 
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